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The limit Gaussian distribution of multivariate weighted func- 
tional of nonlinear transformations of Gaussian stationary processes, 
having multiple singular spectra, is derived, under very general con- 
ditions on the weight function. This paper is motivated by its poten- 
^ C| tial applications in nonlinear regression, and asymptotic inference on 

nonlinear functionals of Gaussian stationary processes with singular 
spectra. 
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1. Introduction. During the last thirty years, a number of papers have 
been devoted to limit theorems for nonlinear transformations of Gaussian 
00 ■ processes and random fields. The pioneering results are those of Taqqu [25, 

26] and Dobrushin and Major [6], for convergence to Gaussian and non- 
^-j- | Gaussian distributions, under long-range dependence, in terms of Hermite 

expansions, as well as Breuer and Major [4], Ivanov and Leonenko [12], 
O ■ Chambers and Slud [5], on convergence to the Gaussian distribution by us- 

ing diagram formulas or graphical methods. This line of research continues 
to be of interest today; see Berman [3] for m-dependent approximation ap- 
proach, Ho and Hsing [9] for martingale approach, Nualart and Peccati [18] 
^ ' (see also Peccati and Tudor [23]) for the application of Malliavin calculus, 
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Nourdin and Peccati [15] in relation to Stein's method and exact Berry- 
Esseen asymptotics for functionals of Gaussian fields, Avram, Leonenko and 
Sakhno [2] for an extension of graphical method for random fields, to name 
only a few papers. The volume of Doukhan, Oppenheim and Taqqu [7] con- 
tains outstanding surveys of the field. Limit theorems for weighted func- 
tionals of stochastic processes, and for processes with seasonalities were 
considered by a number of authors, including Rosenblatt [24], Oppenheim, 
Ould Haye and Viano [19], Haye [20], and their references. Limit theorems for 
nonlinear transformations of vector Gaussian processes have been obtained 
by Ar cones [1]; see also his references. 

In this paper, our main result (Theorem 5.1) states the convergence to 
the Gaussian distribution of the multivariate weighted functionals of non- 
linear transformations ^(^(t)) of Gaussian stationary processes £(£), with 
multiple singularities in their spectra, having covariance function (c.f.) be- 
longing to a parametric family defined in Assumption (A2) below. Here, tp £ 
<^( x ) dx) [see Assumption (A3) in Section 2], where <p(x) = e~ x ^ 2 /\/2tt, 
denotes the standard Gaussian probability density. Specifically, under 
suitable conditions, the convergence to the Gaussian distribution of 

(1.1) Ct = W? 1 f T w(t)^(t)Hdt) 

Jo 

as T — y 00, is obtained for certain ranges of the parameters defining the 
spectral singularities of £; see Assumption (A4) in the next section. For 
each T > 0, 

w(t) = ( Wl (t), . . .,w q (t))', W| = diag(W? T )? =1 , 

(1-2) 

Wf T = f wj{t)v{dt), i=l,...,q, 
Jo 

where, to ensure a finite limit variance, the weak convergence of the family 
of matrix measures associated with w over the intervals {[0,T],T > 0} is 
also assumed, jointly with some restrictions on the boundedness of their 
components in some neighborhoods of the spectral singularities of £; see 
Assumptions (Bl) and (B2) in Section 4. The convergence to the Gaussian 
distribution also requires some conditions to be assumed on the norms of 
components of function w; see condition (B3) in Section 5. 

As commented, the spectral density (s.d.) / of £(t) is assumed to display 
several singularities denoted as H no i sc = {±xo, . . . , ±x r }, with < kq < k\ < 
■ ■ ■ < jCp. In the case where the weak-sense limit of the measures associated 
with the multivariate weight function w is an atomic measure, it is also 
assumed that its atoms S regr = {Si, . . . , 5 n } do not intersect with the singu- 
larities of / (i.e., Si 7^ ±Xj,i = 1, ... ,n,j = 0,1, ... ,r). The convergence to 
the Gaussian distribution then holds with standard normalization. 
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The nature of the limit results obtained depends on the intersection of the 
two spectral point sets H no i se and S reg r- In the discrete case, this phenomenon 
was discussed by Yajima [27, 28] in some other regression scheme. Otherwise, 
different normalizing factors must be derived, and new limiting distributions 
are obtained, for Hermite rank m > 2. Note that the classical noncentral 
limit theorems (Taqqu [26], and Dobrushin and Major [6]) can be viewed 
as particular cases of the general setting considered here, when there is an 
unique singular point in the spectrum of £, with hq = 0, and w(t) = 1. In this 
case, the measure sequence, constructed from the weight function w, is given 
in terms of the Fejer kernel, which tends to the delta-measure with atom at 
zero. Some limiting distributions for the case when the two spectral point 
sets S n oise and H regr are in fact overlapped, in discrete time, can be derived 
from the papers by Rosenblatt [24], Arcones [1], Oppenheim, Ould Haye and 
Viano [19] and Haye [20]. In continuous time, the limiting distributions for 
nonempty set, H no ; se n H regr , can be obtained from the paper of Ivanov and 
Leonenko [13], and the book by Ivanov and Leonenko [12]. This subject will 
be considered in subsequent papers. 

In the derivation of the main result of this paper, Peccati and Tudor's 
central limit theorem [23] (see also Nualart and Peccati [18]), for a family 
of vectors of random variables (r.v.'s) belonging to fixed Wiener chaoses, is 
applied. The outline of the paper is the following. Motivating examples, as 
well as preliminary identities, and conditions needed in the derivation of the 
subsequent results are provided in Section 2. The zero- mean Gaussian ran- 
dom field family considered is embedded into an isonormal process family 
in Section 3. The conditions needed for the weak-convergence (in particular, 
to an atomic measure) of the matrix-valued measures associated with the 
class of vectorial weight functions studied are established in Section 4. The 
asymptotic normality of the corresponding weighted functionals of nonlin- 
ear transformations of zero-mean Gaussian stationary random processes is 
obtained in Section 5. Section 6 provides the final comments, and our main 
conjecture on the work is developed. 

2. Stationary processes with singular spectra. Let us consider simulta- 
neously discrete and continuous time cases in the following development. 
Specifically, for a stationary process £ defined on a complete probability 
space (Q,¥,P), the following notation will be followed: 

((()=((w,i):fix§^l, 

where S = Z, for discrete time t G Z, and S = R, for continuous time t € R. 
Such a process is assumed to be measurable and mean-square continuous in 
the case of continuous time [see also Assumption (Al) below]. 

In the definition of integrals, v(dt) will represent a counting measure in 
the case of discrete time [i.e., v{{t}) = 1, t S Z], and the Lebesgue measure 
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dt, in the case of continuous time [i.e., v(dt) = dt if t E R]. According to this 
notation, the integral 

F 9 (t)£(t)v(dt) 

J 

represents the sum X)fr=i £Ws(*)> f° r discrete time, and the Lebesgue inte- 
gral J g(t)£(t) dt, for continuous time, where g{t) is a nonrandom (measur- 
able for continuous time) function. 

Consider now the following motivating example. 

Example. Let x be defined in terms of the nonlinear regression model 

(2.1) x(t)=g(t,9) + iP(C(t)), t£§+, 

where §+ = R+, for continuous time and S + = N, for discrete time, and with 
g{t, 9) : S + x — > R being a continuously differentiable function of an un- 
known parameter 6 E G C R 9 , consider gi{t,9) = (d/d9i)g(t,6), i = 1, . . . ,q, 
such that 

(2.2) d\ T = \ [ gi (t,9)] 2 u{dt) <oo, r>0,t = l,...,g, 

J o 

and ip(£(t)) represents the noise, with Eip(£(t)) = 0. The least squares es- 
timate (LSE) 0t of an unknown parameter 9 E 0, obtained from the ob- 
servations xit), t E [0, T], or t = 1, . . . ,T, is any r.v. 6? E C , having the 
property 

i-T 

Qt0t)= inf Qt(t), Qt{t)= [x(t) - g(t,r)] 2 u(dt), 

where C is the closure of 0. Let Vg(t,9) = {g\{t,9), . . . ,g q {t,9))' be the col- 
umn vector-gradient of the function g(t,9). We denote dj,(6) = diag(d? T )| =1 , 
where d 2 T , i = 1, . . . ,q, are defined by (2.2). In the theory of statistical esti- 
mation of unknown parameter 9 E C R 9 for the scheme (2.1), the asymp- 
totic behavior, as T — > oo, of the functional 

(2.3) C,T = d- T \9) ! T Vg{t,9)^{t))v{dt), 

Jo 

plays a crucial role, since, under certain number of conditions, the asymp- 
totic distributions of the normalized LSE dx{9){9T — 9), and properly nor- 
malized functional (2.3) coincide, as T — > oo; see Ivanov and Leonenko [12, 13] 

In this setting, an interesting case corresponds to £(t) to be a Gaus- 
sian stationary process with s.d. /(A) displaying singularities at the points 
^noise = {±Xj, j = 0, 1, 2, ... , r}; see (2.5) below. The nonlinear functions 

g(t, 9) = iP cos(it? + 4>), (3 > 0, E R, cp E (— tt, tt],9= (/3, 0, 0) 
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are of particular interest in applications because they themselves also involve 
various seasonalities. 

Let us consider {£(t),t G §>} to be a stochastic process satisfying the fol- 
lowing assumptions: 

(Al) Process £ is a real stationary mean-square continuous Gaussian pro- 
cess with E£(t) = 0, E£ 2 (t) = 1. 
(A2) The c.f. of £ is of the form 

r 

(2.4) B(t) = E[£(0)£(t)] = Y t A j B aj „(t), teS,r>0, 

3=0 

where, for j = 0, . . . , r, 

. . cos(x,-t) 

= (l + t2)aj/2 ' < X < Xi < • • • < X r ,0 < < l,t G S, 

r 

^^• = 1, Aj>0,j=0,...,r. 

3=0 

The c.f. B(t), t G § admits the following spectral decomposition: 

B(t)= [ e iXt f(X)d\, teR, 

where the set A = (— 7r,7r], in the discrete case (t G Z), and A = R, in the 
continuous case (tGM), and the s.d. / in the continuous time is of the form 

r 

(2-5) /(A) = ^i/cWA), AGE, 

3=0 

where, for j = 0, . . . , r, and A G R, 

/ Qj ^(A) = ^[A- (Qj „ 1)/2 (|A + x,|)|A + x j |^- 1 )/ 2 

+K {a ._ l)/2 {\\-K j \)\\->c j \^- i y\ 

with 

2 (l-«)/2 
%JitY{a/2) 

and 

K u (z) = ^J^ s^expj-Vs + i^jds, z>0,^l, 

being the modified Bessel function of the third kind of order v or McDon- 
ald's function. We also note that K_ u (z) = K u (z), and for z I K u (z) ~ 
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Thus, as A — > ±Xj, for j = 0, . . . , r, 

foij,Xj (a) = ^[|A + - h 3 i\\ + *,-|)) 

(2-6) 

+ | A _ x .|^-i(l_/ ii (|A-^|))], 

where 

c 2 (a) 



2r(a) cos(avr/2)' 



^•(|A|) 



r((a j + l)/2) 



r((3-a i )/2) 



A 



l-a 



> r((a 3 - + l)/2) 
4T((3 + a,)/2) 



2 



A— »>0,j = 0,. 



Therefore, the s.d. / has 2r + 2 different singular points [see condition (A2)], 
when kq ^ 0. 

A model with discrete time which satisfies condition (A2) can be obtained 
by using discretization procedure and the formula for s.d. of stationary pro- 
cesses with discrete time of the form 



£ /(A + 2fcvr). 



We will use the same notation for the s.d. in both cases corresponding to 
discrete and continuous time. 

Similar results can be obtained for c.f.'s of the form 

„ . . cosfx,i) m 

R *i,«M = (1 + | f | Pj)tV *j G M >° < <*m < ^ °>J = o, . • . ,r 

(see again Ivanov and Leonenko [13] for details). 

It is well known that the Hermite polynomials H^(x) = (—l) k e x2 ^ 2 -^e~ x2 ^ 2 , 
A; = 0,1,... constitute a complete orthogonal system in the Hilbert space 
L2QSL, y>(x) dx) of square integrable functions with respect to the standard 
Gaussian density if. 

(A3) Assume that the function ip G L 2 (R, <p(x) dx), that is, Eip 2 (((0)) < 
00, and CoO/O = (0)) = 0. 

Definition 2.1. A function -ip £ L 2 (R, tp(x) dx) has Hermite rank 
iJrank(^) = m if either Ci(ip) 7^ and m = 1, or for some m > 2, C\{ip) = 
••• = C7 Tn _i(^)=0,C ro (^)^0. 

(A4) Either (i) i^rank(^) = 1, a > 1/2; or (ii) i^rank(^) = m, am > 1, 
where a = min./ = o,...,r "j, with ay, j = 0, . . . , r, introduced in (A2). 
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Under condition (A3), function ip(x) of Hvaiik.{tp) = m can be expanded 
into a Hermite series in the Hilbert space ^(M, y>{x) dx) 

(2.7) j,(x) = f; ^P-H k (x), 

k=m 

or the process if>(£(t)) admits a Hermite series expansion in the Hilbert space 

L 2 (n,¥,p) 

(2-8) Mffl) = E ^p^(eW), 

where 



Ck{ip) = / ip(x)Hk(x)(p(x) dx, k>0. 
Jr 

3. Some elements of the theory of isonormal processes. In this section, 
we introduce basic notation, elements and results in relation to Gaussian 
Hilbert spaces, isonormal processes and chaos expansions needed for our 
purposes; see Nualart [17]; Janson [14]; Nualart and Peccati [18]; Peccati 
and Tudor [23]; Peccati [21]; Nourdin, Peccati and Reveillac [16], among 
others. 

Definition 3.1. Let H be a real separable Hilbert space. The set of 
r.v.'s X = {X(h) : h 6 H} is said to be an isonormal process on H if X is a 
centered -ff-indexed Gaussian family defined on a probability space (O, F, P), 
and it satisfies 

E[X(h)X(g)] = (h,g) H , h,g£H. 

Let us now consider a real-valued centered Gaussian process £ indexed 
over S = M. By E denote the collection of all finite linear combinations of 
indicator functions of the type l(-oo,t] j with tsl. To embed a real- valued 
centered Gaussian process £ indexed by M into some isonormal process X, 
we introduce a separable Hilbert space H defined as the closure of £ with 
respect to the scalar product 

(3.1) (f,h) H :=J2^c 1 Eas i )C(t j ) 

y 

for given functions / = X)j a i\-oc, Si ] and h = Ylj c j\-oc,tj] m H. Thus, for 
any function h = ^ <k\-oo,u] e ^ ' define 

(3.2) X(/0=J>£(*0- 

i 

Additionally, for any function h£ H, X{h) can be defined as the limit in 
L2(fi,F,P) of X{h n ) for any sequence {/i n } C <S convergent to h in ii". This 
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sequence may be not unique, but the definition of X(h) does not depend 
on the choice of the sequence {h n }. From this construction, process X is an 
isonormal process over H defined as ^(l(_oo,t]) = 

When § = Z, a similar development in terms of sequences leads to the 
definition of an isonormal process from a Gaussian process £ on Z. Now, £ 
denotes the set of all real- valued sequences h = {hi : I £ Z} such that h\ 7^ 
only for a finite number of integers I. The real separable Hilbert space H is 
then introduced as the closure of the set £ with respect to the scalar product 

(f,h) H -.= J2fkhiE{;(k)t;(i) 

k,l 

for given sequences / = {/& -k £ Z} and h = {hi : I E Z}. If h E H, then the se- 
ries ^ g z^z£(0 converges in L2(^,F,P). Thus the centered Gaussian family 
{X(h):h€H}, with 

(3.3) = J>£(0 

zez 

is an isonormal process over iJ. 

Let X be an isonormal process defined on H as before, that is, from a 
centered Gaussian random process £. Let us write Hq(X) =M 1 , and %\(X) 
the closed linear subspace of the set of r.v.'s {A^(/i) : /i £ i/} in the Hilbert 
space L 2 (0,F,P). Thus 

X:H —>Hx(X), 
h^X(h). 

For any n > 2, by T~L n (X), the nth Wiener chaos of process X is denoted, 
that is, the closed subspace of £2(0, F,P) generated by the r.v.'s H n (Y), 
where Y € H\(X), and i?[Y~ 2 ] = 1, with H n denoting, as before, the nth 
Hermite polynomial. Let us now consider the isometry 

(3.4) lX:H° n ^U n {X), 

between the symmetric tensor product H &n , equipped with the norm \/nT|| • 
||#®n, and the nth Wiener chaos Tl n (X) of X. For any h E }{® n ^ I^(h) is 
then defined as I* "(h) := I^(h), with h denoting the symmetrization of h. 
For any g £ and h £ #® n , 

Wmi^In (h)] = 5 mn m\(g,ti) H ® m . 

The pth contraction of g = g\® ■ ■ ■ ®g k £ J?®* and /i = /ii ® • • • <8> h k £ i?' X)fc , 
designated as (7 ® p h, is the element of i^® 2 ( fc_ p) given by 

(3.5) g® p h = (hi,gi) H ■ ■ ■ (h p ,g p )H9p+i ® • • • <8> Vh ®---®h k . 

The definition can be extended by linearity to any element of 
Finally, any r.v. F £ £2(^1 G,P), with <r- field G generated by the r.v.'s 
{X(h),h £ H}, admits an unique chaos decomposition F = X^fcLo^ife (^fc)> 
where /i^ £ if 0fc . 
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From the constructions (3.2) and (3.3) of an isonormal process X from 
a Gaussian process £, respectively, denned over continuous and discrete 
time, H n (X),n > 1, coincides with the nth Wiener chaos associated with 
£, ~H n (,Oi n — 1- Since, by definitions (3.2) and (3.3), %\{X) = T-Li(£), and, 
as stated before, the nth Wiener chaos of process X is the closed subspace 
of £2(0, F,P) generated from the evaluation of nth Hermite polynomial H n 
over the r.v.'s of the space Hi(X) = %i(£). 

The next statement is a convenient, for our purposes, modification of 
Theorem 1 of Peccati and Tudor [23]; see also Nualart and Peccati [18] (in 
the above papers all statements are formulated for positive integers T G 
{1,2,...}, but it is easy to see that one can formulate similar results for 
continuous T > as well) . 

Proposition 3.1. Let {£(t),t G §} be a centered Gaussian process, and 
X is the isonormal process constructed from it as given in (3.2) and (3.3). 
Consider the natural numbers: 1 < n\ < n-i < ■ ■ ■ < < oo, d>2, and the set 
of r.v. 's -K T , nj (£) G H nj (£), where, for T > 0, n T , nj (0 = ^ (fj,r), for certain 
fj s T G H &n i , j = 1, . . . , d, such that 

(3.6) lim ETrl n .{0 = J im n j [ \\fj,T\\'L®n j =1, j = l,...,d. 

Then the following conditions are equivalent: 

(i) For each j = 1, . . . , d, 

lim \\f jtT <g)p f j:T \\ ®2(tw-,,) = 

T— >oo 

/or every p = 1, . . . , n,- — 1. 

(ii) For every j = 1 , . . . , d, 



lim £7[(I^(/ ilT )r] = 3. 

(iii) As T—too, the vector (I,^(/i,t)> • • • j-^(/<2,t)) converges in distri- 
bution to a d-dimensional standard Gaussian vector Nd(0, Id). 

The proof follows from Peccati and Tudor [23], and Nualart and Pec- 
cati [18], considering the fact that H n (£) = ~H- n (X), for any n > 1, with X 
being the isonormal process constructed from identity (3.2), in the continu- 
ous time case, and, similarly, in the discrete time case, from equation (3.3). 

Corollary 3.1. Assume that conditions (3.6) and (i) or (ii) of Propo- 
sition 3.1 are satisfied for r.v. 's 

(3-7) 7r T , ni (£)= f r Tij (t)H nj (at)Hdt), 

Jo 
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where, in the case of continuous time, it is also assumed that rTj(t) E 
C([0, oo)), for T > 0, and j = 1, . . . , d. Then, the vector 

(3-8) tt T4 (0 = (^rT,!^)^^^))^^),...,^^^)^^^))^^) 

converges in distributions, as T — )■ oo, to a standard Gaussian vector tt<i ~ 
N(0,I d ). 

Proof. In the case of continuous time, since = -X"(l(_oo,t])> 

where denotes the isometry introduced in (3.4). Therefore, for rTj(t) E 
C([0,oo)), T>0 and for j = 1, . . . ,d, 



*T, nj (£) = £ r Ttj (t)H nj (£(*)) dt = I* r T>j (t)l 



, 3 ,dt). 

(— oo,t\ 



Thus, considering in (iii) of Proposition 3.1 

fj,r(si , . . . , s n .) = J r T ,j (*)lf_ n 4, i4] (si , ■ ■ ■ , s n . ) dt 

for j = 1, . . . , d, we obtain the desired result. 

Similarly, for the case of discrete time, we have, from (3.3), 

(3.9) X(6. tl )=Z(l), lez, 

where, for each I E Z, S.i denotes the Kronecker delta function, that is, 

1, if i = I, 
0, ifi^l,ieZ. 
Therefore, 

H nj m) = H n] (X(6, l )) = lX(8f > ;>). 

Proposition 3.1 (iii) is then applied, considering 

T j 

/ jiT (mi,...,m„.) =^2 r T,j{l)Y[ S ™ l ,i, mi,...,mj eZ 

1=1 i=l 

for j = 1, . . . ,d. □ 

4. Spectral measures of weight functions and admissible spectral den- 
sities. Let us first establish some results on weak-convergence of matrix- 
valued measures, given by 



(a n Uriw w 3 T {X)w l T {X) d\ . 

(4.1) /4(dA)= 1 =, j,l = l,...,q, 



^f A \wUX)\ 2 d\J A \w l T (X^d\ 
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where 

cT 



w 3 T (X)= e ltx Wj {t)v{dt), j = l,...,q, 
Jo 

and the functions Wj(t), j = 1, . . . ,q, are, as before, the functions (1.2) in- 
volved in the definition of the random vector (1.1). 

(Bl) Assume that the weak-convergence fi T =^> /x, when T — > oo holds, 
where n T is defined by (4.1) and /i is a positive definite matrix measure. 

The above condition means that an element fjp 1 of the matrix-valued mea- 
sure n is a signed measure of bounded variation, and the matrix n(A) is 
positive definite for any set A€ A, with A denoting the a-algebra of mea- 
surable subsets of R; see, for example, Ibragimov and Rozanov [10]. 

The following definition can be found in Grenander and Rosenblatt [8], 
Ibragimov and Rozanov [10] and Ivanov and Leonenko [12]. 

Definition 4.1. The nondegenerate matrix- valued measure fi(d\) = 
{//^(dA)}* l=1 is said to be the spectral measure of function w(i). 

Definition 4.2 (Ibragimov and Rozanov [10]). The s.d. / is said to be 
/i-admissible if it is integrable, that is, all elements of the matrix 

/ /(A)M(dA) 

J A 

are finite, and 

(4.2) lim / f(\)ti T (d\) = [ /(A)fi(dA). 

Let us introduce two conditions on the s.d. / that guarantee its /i- 
admissibility. These assumptions are related to basic conditions on the c.f. 
and s.d. (A2). In the following, J denotes one of the three sets: 

{-r, . ... -1, 0, 1, ... ,r}; {-r, . . . , -1, 1, . . . , r}; {0}. 

We formulate the following condition for a set J = {—r, . . . , — 1, 0, 1, . . . , r}. 
(I) The s.d. feC(A\ {xj,j e J}), with 

x_j = — Xj, j = 0, 1, . . . ,r,0 < xo < >c\ < ■ ■ ■ < x r 

and, for j = 0, 1, . . . , r, 

lim /(A)|A-x j | 1 -^ = aj >0, 

A— ¥>cj 

(Xj G (0, 1), j G J; a-j = aj,a^j = aj. 
We obtain from (4.3) that, for any e > 0, and j G J, there exists 8j = Sj(e), 

a,- + e 



such that for | A — x,- 1 < 5j 



/(A)< 



x/-«, 
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Then, for |A — Hj\ < 5j, we have the following: 



{\: f(\) > c} CVj{c) = < \:\\ - xj\ < 



Moreover c must satisfy the inequality 



aj +e 



and equivalently, 

c>^±f-= Cj (s). 

S j 

(II) Let £o > be fixed. There exists Co = maxj e jCj(eo) 5 such that for 
c > c , 

(4.5) {\:f(\)>c}c\JV j (c), 

where Cj(e) are defined by (4.4). 

It is easy to see that for sufficiently large c (say, c > Co), the neighborhoods 
Vj(c), j £ J, in (4.5), are nonoverlapping, and 

|^( c )|;o 

as c — > oo. 

Note that the function (2.5) satisfies conditions (I) and (II). 
(B2) For T sufficiently large (say, T > T ), 

(4.6) Wff max |it>r(A)| < % < oo, j 6 J, i = 1, . . . , q. 

X€Vj{co) 

In condition (B2), one can assume that (4.6) holds only for j = 0, 1, . . . ,r, 
since V-j(co) = -Vj(c ), j = 0, 1, . . . , r. 

Theorem 4.1. Assume that conditions (Bl), (B2), as well as (I), (II) 
are satisfied, and the s.d. f is integrable with respect to the spectral measure 
li, then the s.d. f is fi-admissible. 

Proof. For c > Co, we consider 

/ c (A) = /(A)l (/(A)<c) (A)+cl (/(A) > c) (A). 
Then, for k, I = 1, . . . , q, 

/(A)4''(dA) - / f(X)^ l (dX) 
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< 



/(A)4''(dA)- / f c (\W(d\) 



fc/ \\ ,,k,li 



(4.7) 



+ 



+ 



n\)^\d\) - 1 n\)^ i (d\) 



n\)^ i {d\)- / f(\)^ i (d\) 



= I k 1 \T 1 c) + I k 2 \T 1 c) + I k /{c). 
By Assumption (Bl), for any complex numbers z = (z\, . . . , z q ), the function 

g 

M Z (A) = Yl ^\ A ) z ^l ^ 0, A G A, 

k,l=l 

is a measure. Thus, by Lebesgue's monotone convergence theorem, 

(4.8) / f(\)M z (d\) — > / f(X)M z (dX). 

Note that the diagonal elements fi k ' k and are measures; thus if only Zk 
and zi are nonzero among 2; = (zi, . . . , z q ), we obtain from (4.8) that 

(4.9) [ (f(X)-f c (X))( f i k ' l (dX)z k z l + ^ k (dX)z l z k )^0, c^oo. 
J A 

Note that fj ,k = fj, k > 1 , and choosing, for instance, z k = z\ = 1, we have from 
(4.9), 

' (/(A) - / C (A)) Re(fi k ' l )(dX) — > 0, c -> 00. 

L 

If we choose z& = 1, #j = — i, then 

(/(A)-/ c (A))Im(/^)(dA)^0, c^oo. 



Thus 



lim I kJ (c) =0. 



For a fixed c, we obtain, from condition (Bl), that 



lim / 2 M (T,c) = 0. 

On the other hand, under the conditions assumed in this theorem, for 
T>T , 



{A:/(A)>c} 
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(4.10) 

<f V^,, / /(A)dA-X), 



when c — > oo. Thus, for any e > and T > Tq, one can choose c\ = ci(e) > Co, 
such that for c > ci, we have I^ ,l (T,c) < e/3. Then, once can take C2 = 
C2(e) > co, such that for c > C2, we have I*' (c) < e/3. 

Let us now fix c = max(ci,C2); then, there exists T\ = T\(e) > To, such 

□ 



that for T > Tq, 1%' (T, c) < e/3, and the left-hand side of (4.7) is less than e 



5. Central limit theorem for weighted functionals. This section provides 
the asymptotic normality as T — > oo of the vector (1.1), that is, we will prove 
that the vector £ T converges in distribution (=$■) to some Gaussian vector £. 
Thus, for any z £ M 9 , we prove that z) ==> (£, z), as T-> oo. Denoting, 
for z = (zi, . . .,z q ), 

1 

Y^ZiW^Wiit) = R T (t, z) = R T (t), 

i=l 

from (2.8), we have 

(C T ,z}= [ T mt ))R T (t)u(dt) = f2^1 j T R T {t)Hj{^(t))u(dt). 
Jo J- Jo 

In the derivation of the proof of our main result, the following additional 
conditions are required: 
(B3) For T>T , 

(5.1) W7^su V \w i (t)\<k i T- 1 l' 2 1 i = l,...,q, 

where the supremum is taken over t in the interval [0, T], in the case of 
continuous time, and over t in the set {1,...,T}, in the case of discrete 
time. 

Let / ( * 1} (A) = /(A), and for j > 2, 

f j 
/* 0) (A)= / /(A-A 2 \j)T\f(\i)d\ 2 ---d\j, 



the jth convolution of the s.d. /(A). 
(C) The matrix integrals 



/ /*«)(A)/x(rfA), j>l, 
Ja 



are positive definite. 

We now proceed the formulation of our main result. 
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Theorem 5.1. Suppose that conditions (Al)-(A4), (B1)-(B3) and (C) 
are fulfilled. Then, the r.v. C, T in (1.1) converges in distribution, as T — > oo, 
to the Gaussian r.v. £ with zero mean and covariance matrix 

(5.2) S = 2vr ^ -J^l / f^(XMd\), 

where /x is the weak-sense limit of the family of matrix-valued measures 
introduced in (4-1) and associated with the weight function w(t) in (1.2), 
given from functional (1.1). 



In the proof of the above theorem, the following identities will be applied 
jointly with Lemma 5.1 formulated below. Specifically, from the orthogonal- 
ity of Hermite polynomials, we obtain 



E(C T ,z) 2 =J2 



C;(V) 



1=m 

(5.3) 

00 C|(V) rT rT 



— (J [IP) f f 

y^_iw / R T (t)R T ( s )Bi{t-s)v{dt)v{ds) 



We will prove the asymptotic normality of (1.1) under condition (A4)(i). 
The proof under condition (A4)(ii) is even simpler. 

By conditions (A2) and (A4)(i), for j > 2, all the convolutions /*W are 
bounded and continuous functions, and by (Bl), 



<>tU,z) = j\ l T C B\t-s)R T {t)R T {s)v{dt)v{ds) 
Jo Jo 



= 2^j\Jj^\x)i^^\X)z k z}jdX 
— > 2ttj! / f<i\X)m z (dX) = a\j,z), j>2, 

T->oo J A 

where m z {dX) = YX,i=\l jk,l ( d)K ) z k z i- 

Under condition (A2), from Theorem 4.1, we obtain for j = 1, 

(5.5) lim <4(1,z) = 2tt / f(X)mJdX) = a 2 {l, z). 

T^oo J A 



(5.4) 
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-i 2 



a 2 (j,z) = a 2 (z). 



lim £<C r ,z) 2 = ]T 
In Lemma 5.1 below, we will consider the following decomposition: 

TT = (Ct, ^) = 7^(d) + Tj<(d) 

(5.7) 



E+ E 

<3=1 j=d+l/ 



cm 

f- Jo 



RT(t)Hj^(t))u(dt). 



Lemma 5.1. Suppose that conditions (A1)-(A4) and (B1)-(B3) hold. If 
for any d>l, as T — > 00, Tr(<f) =^ ~ iV(0, cr^(z)), where 

d 



(5.8) 



°3 



(*)=£ 



C,-(^) 



then t t ^-t~N(0,o- 2 (z)). 



Proof. Note that E[(Tfp(d)) 2 ] ->• 0,d — >■ 00, uniformly in T. Really, by 
condition (B3), 



(5.9) \R T {t)\ 



t=l 



<T _1/2 ||z||||fc||, k = {k u ...,k q ). 



Then, under (A4)(i), as d— >oo, 

c 2 W 



_.~T7, J- Jo Jo 



j=d+i 



<r- 1 lUll 2 ||fc|| 2 



B 3 {t - s)R T (t)R T (s)v(dt)v(ds 

[ [ B 2 (t- s)u{dt)v{ds) V 
Jo Jo ■_, , 



c 2 M 



j=d+l 



„ 00 

< ||z|| 2 ||fc|| 2 / B 2 {t)v{dt) 

J R j=d+l 



c?(V) 



(3(d) — >0, 



since by Parseval's identity, 

00 C 2 (ib) 



Thus, for any e > 0, uniformly in T, 

p{|4(d)|> £ }<^ 



0, d— >oo. 
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For any e > 0, and d > 1, we then obtain 

(5.10) lh^P{T T <x}<<S> d (x + e) + ^, 

where is the distribution function of a Gaussian r.v. with zero mean and 
variance cr^(z). 

Also, for any e > 0, and d > 1, as T — > oo, the following inequality holds: 

(5.11) ljmP{T T <x}><f> d {x-e)-^. 

e 2 

If d— > oo, we obtain, from equations (5.10) and (5.11) that, as T — > oo, 

$oc(x - e) <hmP{T T < x} < limP{rr < x} < &oo(x + s), 

where $oo is the distribution function of a Gaussian r.v. with zero mean 
and the variance a 2 (z) given by (5.6). Thus, if e — > 0, limT-^oo P{tt < x] = 

$oo(4ieK. □ 

Now we are in position to derive the proof of Theorem 5.1. In such a 
proof, we will check condition (i) of Proposition 3.1, but the proof can also 
be developed from the verification of condition (ii) in Proposition 3.1, using 
diagram formula. We place this proof into Appendix, due to its method- 
ological interest in relation to the approach it presents for the analysis of 
nonregular diagrams, providing the classification of their levels into recipi- 
ents and donors. 



Proof of Theorem 5.1. From Lemma 5.1, it is sufficient to show the 
asymptotic normality of the r.v.'s T-r(d). Consider then the r.v.'s 



(5.12) 7T T4 (C) 

where 
(5.13) 



rr,i(t)i2i(£(t)Mcft), 
Rr{t) 



T \ l 

r T Jt)H d (£(t))v(dt) 



o-(j,z) 



j = l,...,d. 



The proof will follow from the application of Corollary 3.1, after checking 
condition (i) of Proposition 3.1 for the random vector 7TT,d(£,) defined by 
(5.12) and (5.13). From Theorem 1 and equation (5.4), 



E 



(5.14) 



<4(j» 

o 2 {j,z) 



1. 



T-^oo,j = l,...,d. 
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Now, 7Pr,d(£) ==> iTd ~ A/"(0, Id), T — > oo, if and only if 
lim \\fj,T<S>pfj,T\\ H ®2(j-p) =0 

T— yoo 

for p = 1, . . . ,j — 1, 2 < j <d, where 

pT 3 

f jtT (si,...Sj)= / R T (t)Y[l(-oo,t](si)dt. 
•'° 1=1 

We first check the convergence to zero of contractions in the continuous 
time case. The pth contraction is computed by applying formula (3.5) with 
k = j as follows: 

fj,T ®p fj,T(xi, . . .,X 2 j-2p) 

= [ [ R T (t)R T (s)B(t-s)x---xB(t-s) 
Jo Jo P 

j j 

(5.15) x Y[ \-oo,t]{xi) II l(-oc,s]{xi)dsdt 

i=p+l l=p+l 



j j R T (t)R T {s)B p (t-s) 
Jo Jo 



j-p "2j-2p 

x n^- 00 .*]^) n k-oo, S ](xi) ds dt. 
i=i i—j— p+i 

The norm of the pth contraction (5.15) in the space H® 2 ^~ p ^ is then given 

by 

\\fj,T <8>p fj,T\\H®*ti-p) 

rj~\ rj^-i r^i 

= [ [ [ [ J R r (t 1 ) J R T ( Sl )^ T (t 2 ) J R T (s 2 ) 
Jo Jo Jo Jo 

x Bj- p (h - t 2 )B' J -P{ Sl - s 2 )B p (t l - Sl ) 
x B p (t 2 - S2) dsi ds2 dti dt2- 
By condition (B3), for 2 < j <d and p = 1, . . . ,j — 1, 

\\fj,T ®p fj,T\\ 2 H m(i-v) 



< 



< 



| Z ||4||M|4 f-T r T f-T r T 



rp2 



J<2 



[ [ [ [ \BP(t 1 -8 1 )&>(t 2 -8 2 )\ 

Jo Jo Jo Jo 

x \B j ~ p (ti - t 2 )B j - p ( Sl - S2)\ds 1 ds 2 dt 1 dt2 

r r r [ r \B(t 1 -s 1 )B(t2- S 2)\ 

Jo Jo Jo Jo 



HI 4 !!* 1 " 4 rT rT rT rT 
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(5.16) x |J3(tx -t 2 )B{s 1 - s 2 )\dt 1 dt 2 dsids 2 

n ||4||M|4 r T r T r T -\ r T 



< 



J>2 



I I III [B 2 (h-s 1 ) + B 2 (t 1 -t 2 )]dt 1 
Jo Jo Jo 1 Jo 



x \B(t 2 - s 2 )B(si - s 2 )\ ds\ ds 2 dt 2 



<4 z 4 k 



4111.114 







B 2 (t 1 )dt 1 / \B(t 2 )\dh 



o 



xT 2 f [ \B( Sl -s 2 )\d Sl ds 2 . 
Jo Jo 

In (5.16), as T — > oo, 



(5.17) 



f\B(t 
Jo 



2 )\dt 2 = 0(T 



l-a\ 



T~ 2 f T [ T \B(s 1 -s 2 )\ds 1 ds 2 = 0(T- a ). 
Jo Jo 



From condition (A4), in the case considered of Hermite rank m = 1, we 
have a > 1/2, and therefore, from (5.16) and (5.17), we obtain for j > 2, 



(5.18) 



lim \\fj,T® P fj,T\\H®2U-p) =0- 

T— »oo 



The proof in the discrete time case can be similarly derived in terms 
of definition (3.3) of isonormal process X, considering the counting measure 
is(-)- Specifically, from (3.9), for T > and 2 < j < d we consider the sequence 
of kernels 

T j 

fj,T{mi, . . . ,rrij) = ^Rrjl) ]~[ S mii i, mi, . . . , rrij € Z. 

1=1 i=l 

For p = l,...,j— 1, the pth. self-contraction of this kernel is given by 
fj,T ® P fj,T(mi, . . . ,m 2 j_2 P ) 

T T j-p 2j-2p 

q=l 1=1 i=l i=j—p+l 

Therefore, since 

\\fj,T ®p fj,T\\ H ®2U- P ) 

= E E E E flr(«Wi)^(*)*r(0 

5=1 fc=l Z=l i=l 

x ^'- p (g - k)B j ~ p {l - i)B p (q - l)B p (k - i), 
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in a similar way to the continuous time case, we obtain 

\\fj,T ® v fj,T\\ 2 H ®2(j-v) 



(5.19) 



y2 



UII 4 IIA; 



q=l k=l 1=1 i=l 



x \Bi- p (q-k)B j - p (l-i) 



q=l k=l 1=1 i=l 



< 4IU 



|4||L||4 



-3=1 
T T 

r-EEi^ 



5>(*oi 



,fc=i 



Z=l i=l 



Thus, as T — >• oo, 



Ei^)i = o(^- 



(5.20) 



fe=i 



T T 



r^EEl^-OHocr-"). 



/=1 i=l 



Again, from condition (A4), and equations (5.19) and (5.20), Proposi- 
tion 3.1(i) holds, and the convergence to the Gaussian distribution follows. 



Example (Continuation). Consider now model (2.1) with nonlinear re- 
gression function 



N 



(5.21) 



9(t,0) = ^2(A k cos ip k t + B k sin tp k t), 



k=l 



where = (A 1 ,B 1 ,<p 1 ,...,A N ,B N ,<p N ),C% = Al + Bl>0,k = l,...,N ) 0< 
(pi < ■ ■ ■ < ifN < oo. In this case, q = 3N, function g(t, 9) then has a block- 
diagonal measure n(d\) (see, e.g., Ivanov [11]) with blocks 



Kfc ipfc fa 
-i-Pk, K k 7 fc , 

fa 7k, K k , 



k = l,...,N, 
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where 

ZC k 

Here, the measure = K^dX) and the signed measure p k = Pk(dX) are 
located at the points ±<pk, and /%({±^fc}) = i,pfc({±^}) = ±i. We then 
have 

a d 

93k-2{t,9) = ——g(t,9) = cos ip k t, 93k-i{t,0) = -^-gfaO) = shnp k t, 
oA k dB k 

d 

93k{t,0) = - — g(t,0) = -A k tsmcp k t + B k tcostp k t, k = l,...,N. 
o<Pk 

It is easy to see that if the s.d. / satisfies (I), and jtj ^ (f k ,j = 0, 1, . . . ,r, 
k = 1, . . . , N, one can find a neighborhood Vj{c$) of the point x,-, for j = 
0,1,..., r, which does not contain the points ip k , k = 1, . . . , N. Thus, for 
T > Tq, the following condition holds: 

W^ 1 max \w l T (X)\ < k^T" 1 / 2 , j G J; i = 3fc-2, 3k- 1, 3&; k = 1, . . . ,N. 

Aev,(co) 

In relation to the considered function w(t) = Vg(t, 9), the measure /j, J ^(dX) - 
p,j.(dX,9) approximates, in the weak sense, the spectral measure fj,(dX) = 
{p-* 1 (dX)} Q j l=1 of the nonlinear regression function g(t,9) [see (2.1)], where 



p J T {dX,9) = j,l = l,...,q, 

y/jM(X,9)\ 2 dXj A \g l T (^)\ 2 dX 

g j T {X,9)= [ T ^ tx gj (t,9)v(dt), j = l,...,q 



and gj(t, 9) defines the jth component of w(i) = Vg(i, 9), for j = 1, . . . , q. 

If the s.d. /(A) satisfies condition (II), then / is /x-admissible, and the 
block-diagonal matrix J^f(X)fx(dX) consists of the blocks 
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It is easy to see that, for the function g(t,8) given by (5.21), the matrix 
is a block-diagonal with blocks of the form 

jv/3 \ 
■2C k 



( 



j=m 



J- 



1 







IB ^ 



-A, 



2Cfc 



1 



fc = l,...,]V. 



6. Final comments. This paper addresses the problem of Gaussian limit 
theory of weighted functionals of nonlinear transformations of Gaussian sta- 
tionary random processes £ having multiple singularities in their spectra. 
The general case where the Fourier transform of the weight function also 
displays multiple singularities in the limit, which do not coincide with the 
singularities of the spectral density of £, is also covered here. This subject 
has several applications in asymptotic statistical inference. We are especially 
motivated by its application in the limit theory of nonlinear regression prob- 
lems with regression function and errors having multiple singularities in their 
spectra. This actually constitutes an active research area, due to the exis- 
tence of several open problems and applications. Note that, although here 
we have considered the parameter range 

a = min aj > 1/2, 
j=0,l,...,r 

which, in particular, allows us to consider long-range dependence models. 
Our conjecture is that the Gaussian limit results hold for ay 6 (0,1), j = 
0,1, ... ,r. The proof of this conjecture will lead to a general scenario where 
most of the limit results derived for random fields with singular spectra (see 
Taqqu [25, 26]; Dobrushin and Major [6]; Nualart and Peccati [18]; and the 
references therein) can be obtained as particular cases. 



APPENDIX: PROOF OF THEOREM 5.1 BASED ON DIAGRAM 

FORMULA 

As before, we will prove this result for Hermite rank m = 1. To show 
the asymptotic normality of the r.v.'s tt((1), consider the r.v.'s iTT,d(^) and 
r T ,j(t), j = l,...,d, defined by (5.12) and (5.13). 

We will check condition (ii) of Proposition 3.1. Then, from Corollary 3.1, 
nT,d(O^K d ~N(0,I d ), that is, T T (d) => T d ~ N(0,a 2 d (z)), as T^oo. 

We apply diagram technique for proving condition (ii) of Proposition 3.1. 
Let us first introduce some definitions. 
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A graph T = T(li, ...,l p ) with l\ + ■ ■ ■ + l p vertices is called a diagram of 
order . . . , l p ) if: 

(a) the set of vertices V of the graph T is of the form V = Uj=i Wj > where 
Wj = {(j, l):l<l< lj} is the jth level of the graph T, 1 < j < p (if lj = 0, 
assume Wj = ) ; 

(b) each vertex is of degree 1; 

(c) if ((ji,h), (j2,h)) G r, then j\ ^ j 2 , that is, the edges of the graph Y 
may connect only different levels. 

Let L = L{1\, . . . , lp) be a set of diagrams T of order (li,...,l p ). Denote by 
Zr the set of edges of a graph T £ L. For the edge w = ((ji,h), £ Zr, 

ji < 32, w e set d\{vo) = ji,d 2 (w) = 32- We call a diagram T regular if its levels 
can be split into pairs in such a manner that no edge connects the levels 
belonging to different pairs. We denote by L* the set of regular diagrams 
L* QL(h,...,l p ). Up is odd, then L* = 0. 

The following lemma provides the diagram formula; see Taqqu [26] , Lem- 
ma 3.2 or Doukhan, Oppenheim and Taqqu [7], page 74, or Peccati and 
Taqqu [22]. 

Lemma A.l. Let (£1, . . . ,£ p ),p > 2, be a Gaussian vector with E£j = 0, 
E£j = l, E Zi€j = B (hj),h3 = P, an d let H h (u), . . . , H ip (u) be the Her- 

mite polynomials. Then 

(a.i) 4n>^)l=E n B(d 1 (w),d 2 (m)). 

lj=i J reiraeZp 

From (A.l), we obtain, for p = 4, h = h = h = h = j, T = and 
= (£(*iU(t2),£(t3), £(t 4 )), 

/■T /-T /-T /-T 4 

(a.2) = / / / / n^(*i) 

JO JO JO JO „•_ 1 



1=1 
x £" 

We then have 

(0 



,t=i 



u{dt\)v{dt2)i l {dt^)u{dti). 



T fT -\ 2 

B(ti - t 2 )R T {ti)RT{t2)u{dti)v{dt 2 



ff 4 ! 1 .^ uo jo 

(A.3) 

<7 4 (l,z) 



24 IVANOV, LEONENKO, RUIZ-MEDINA AND SAVICH 

For j > 2, the sum in (A.l) is split into two sums corresponding to regular 
and nonregular diagrams, 

£•••=£•••+ E -. 

reL reL* reL\L* 

and the right-hand side of (A. 2) is split into two these parts, as well. 
Analysis of the regular diagrams: 
We have 

* 

(A.4) £(T) = £ F r (T), 

reL* 

where 

F r (T) 



pT pT pT pT 4 

(A.5) = / / / / Hrrjfc) 

JO J Jo JO 



1=1 

X 



Each regular diagram r € L* consists of 4 levels of cardinality j. There 
are only 3 subdivisions of the 4 levels into pairs, and in each pair the vertices 
can be connected by j\ ways. Thus, there is only 

\L*\ =3(j!) 2 

regular diagrams, and, in this case, sum (A.4) is subdivided into product of 
pairs of integrals 

E( r ) = 4r^rf / T / T ' Bi{t 1 -t 2 )RT(t 1 )RT(t2Mdt 1 )v(dt 2 )) 2 

(A.6) 

Analysis of the nonregular diagrams: 
First, we consider 

(A.7) £(T)= £ F r (T), 

reL\L* 

where i*r is defined as in (A.5). We now prove that lim^oo ^2(T) = 0. Then, 
the assertion of the theorem will follow from (A. 3) and (A.6). 
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From (5.9), 

rj-i rjp rj~i rj-i ^ 

(a.8) >< / / / / n n 

./o ./o ./o ./o i= i w6Zr>dl(w)=i 

x u(dti)v{dt2)y{dt'i)v{dti) . 

Let gr(^) be the number of edges w £ Zr, such that d\{w) = i. Then, for 
ar(») > 1, 

y n i^(*i-*d 2 M)k(^) 

(a.9) E f i^-<4w)r (, M*) 

Jo 



If f/rW = 0, the integrals regarded to these variables (£4, and possibly £3), 
in the left-hand side of (A.9), give a contribution in the form of a multiplier 
of T in the estimate (A.8). 



Definition A.l. The level i of a nonregular diagram r G L \ L* is said 
to be a donor, if f/r(z) > 1 5 an d a strong donor, if <?r(0 = i- The level i of 
a nonregular diagram T £ L\L* is said to be a recipient, if it is not donor, 
that is qr(i) = 0. 



Let p s d be a number of strongly donor levels, and p r bea number of recip- 
ient levels. Obviously, level 1 is a strong donor, while level 4 is a recipient. 
If Psd — 1) then p r = 1, while if p s d = 2, then p r = 2. 

Formulas (A.8) and (A.9) then imply 

II ll 4 lliM|4 JL fT 

(A.10) \F T (T)\ <2 4 -^ l|z '' J 1 r 2 TT / is(£)| 9r(i W£). 

Since j > 2, and a > 1/2, for a strong donor level i with qr{i) = j, 

pT poo 

(A.ll) / \B(t)\ j u{dt)< / [B(t)] 2 u(dt) < 00. 



26 



IVANOV, LEONENKO, RUIZ-MEDINA AND SAVICH 



Thus, for the recipient levels (qr(i) = 0) and the strong donor levels 
(qr(i) = j), we obtain 

(A.12) f 1 \B(t)\ qr{i) v(dt) < CoT 1 -^, 

Jo 

where 

2(i) = ^M C = max(l, J™ B 2 {t)v{dt)^j. 

Let now < qv{i) < j', that is, level i is a donor, but not strong donor, 
and then 



r \B(t)\ qr{i) v{dt) 

Jo 



i 



+ 



B(t)\ qr{i) v{dt) 



rpl-aq T {i) _ 2 

(A.13) + H TV 

1 - aqr{i) 



aq r (i) , T 1 -*"*® _^_ z{ih 



+ i ttt = o(T l 



otqr{i) — 1 1 — «<?r(^) 
since agr(z) = cxjz(i), and aj > 1. We will show that 

4 

/i = 2-£*(i)=0. 



Indeed, 



J>(i) = l + 



gr(2)+9r(3) 



i=l 

and <7r(2) + <7r(3) = J, since |Zr| = 2j. 

Formulas (A.12), (A.13) and (A. 9) together with (A. 10) then imply that 

(A.14) \F r (T)\=0(l), T^oo, 

when p S( j = p r = 2, and 

(A.15) \F T (T)\=o(l), T^oo, 

when < qr(i) < j, for i = 2, 3 (/9 s< j = p r = 1). 

The estimate (A.14) is not exact. Thus, let us consider again the case of 
nonregular diagram T, which has 2 strong donor levels, and the remaining 
2 levels are recipients. The recipient level 3 takes edges from the strong 
donor levels 1 and 2, while level 2 does not supply level 3 in full. Let us 
permutate levels 2 and 3, and denote this permutation by 7r, that is, 7r(2) = 3, 
7r(3) = 2, and, from the level 7r(3) to the level 7r(2), there are less than j 
edges. Moreover, from the level 7r(3) there is no edges down, except the edges 
which connect vr(3) with vr(2), since level vr(3) took all edges from the top, 
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that is, 

<?7rr0r(3)) = q w r{2) <j, 

where 71T is a nonregular diagram, taken from T by permutating the levels 
2 and 3. Note that this permutation does not change the value of integral 
defining Fr(T) in (A. 7), since it is equivalent to the renaming of the variables 
t2 and £3. From (A. 15), we then obtain 

(A.16) \F r {T)\ = \F nT (T)\^0, T^oo. 

The assertion of this theorem then follows from equations (A.14)-(A.16). 
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